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ABSTRACT 



o ! The current modelling of single microlensing light curves neglects the 

r^ I possibihty that only a fraction of the light is due to the lensed star, the 

remaining being due to a close, unresolved blend, which may be related or 
CN ! unrelated to the lens. 

Unfortunately, the effects of blending are significant as all microlensing 
^ I experiments choose very crowded fields as their targets. In this paper we point 

Q\^ I out a strong degeneracy of the fitting procedure which makes it practically 

^ I impossible to detect the presence of a blend by purely photometric means, 

O ! except in a small part of the parameter space. Some blends may be detected by 

Q,. I astrometric means, but the majority have to be corrected for statistically. The 

"^ I luminosity function reaching well below the ground based detection limit (with 

Qh| the HST) would be very helpful. The statistics of binary stars in the target 

O I population is also important and this could be determined with the repeating 

^ I microlensing events. 

If no correction is made then the event time scales, the lens masses, and the 
optical depth are all systematically underestimated. 



Subject headings: galaxy: structure - ISM: extinction - photometry 



1. Introduction 

Almost all model fitting of single microlensing light curves to the data is done with 
a 4-parameter curve which assumes that the stellar image is not blended (cf. Paczynski 
1996, and references therein). However, the first two double lenses: OGLE 7^7 (Udalski et 
al. 1994b) and DUO ^^2 (Alard et al. 1995) were found to have apparent images made of 
at least two objects: the lensed star and another unrelated star. Also, at least one single 
event, OGLE ^5, was found to be strongly blended (Mao 1995, Alard 1996b). According to 
Alard (1996b) most DUO events are likely to be blended, and fitting data with unblended 
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light curves introduces a systematic bias in the estimate of event time scales and optical 
depth. Other effects of blending were considered by Nemiroff (1994), Di Stefano & Mao 
(1995), and by Buchalter & Kamionkowski (1996). 

Another issue is the lensing of very faint stars which become detectable only while 
lensed (Nemiroff 1994). Such objects are missed in searches like DUO (Alard 1996a), 
EROS (Auburg et al. 1993), MACHO (Alcock et al. 1993) and OGLE (Udalski et al. 
1993), but they are detectable in the searches which use the "image subtraction" technique, 
hke AGAPE (Bouquet et al. 1996) and COLUMBIA- VATT (Crotts & Tomaney 1996, 
Tomaney & Crotts 1996, 1997). An elaborate theory of this technique was developed by 
Gould (1996). 

The aim of this paper is to demonstrate some of the practical consequences of blending, 
and the limitations imposed by unknown blends on the determination of lens parameters. 



2. The Model 

Let us consider an idealized situation with a variability due to microlensing superposed 
on arbitrary constant background of whatever nature: the sky, or the crowded field made 
of many stars or nebulae. We adopt an approximation according to which in a given 
small aperture (it may be profiled as the PSF, the Point Spread Function) there is a well 
measured level of brightness long before and long after the microlensing event: 

Fo±AoA^o-^/^ iVo » 1, (1) 

where Aq refers to a standard deviation of a single measurement, and the number of 
measurements Nq is very large. Therefore, Fq as well as Aq are known very accurately. 
However, we do not know what fraction of Fq is due to the star which is microlensed while 
it is in its normal, unlensed condition. 

Let the flux from the lensed star be given as 

Fs = FsoA{t), Fso = fsFo < Fo, (2) 

where the magnification due to microlensing is 

U + Z 2/,\ _ 2 . ft ~ ir 

u (m^ + 4) ' 

where Umin is the impact parameter in units of Einstein ring radius, tmax is the time at 
which maximum magnification is reached, and to is the time it takes the lens to move with 
respect to the source by one Einstein ring radius (Paczyhski 1996, and references therein). 



m = — T^^, u\t) = «L„ + ^— ^ , (3) 
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Given a set of photometric measurements: {Fk,tk), "^ = 1, 2, 3, ..., n, we would 
like to determine four parameters: F^q, Umin, tmax, ^o- In principle the value of the fifth 
parameter: Fq has to be determined as well. However, in practice a good microlensing 
search is conducted over a time interval much longer than the microlensing time scale to, 
so we may simplify our task by adopting Fq as known. In a common fitting procedure it is 
assumed that F^q = Fq, and Fq has to be determined from the same set of measurements as 
the other three parameters: Umin, imax, ^o- 

We assume that the errors (standard deviations) of the measurements scale as 

A.= (^)'^'ao, (^) = (1-/.) + M(4), /.^^, (4) 

where Aq is the error of a single measurement at intensity Fq, Jq is the fractional intensity 
of the lensed star far from the microlensing event, and the magnification A{tk) is given with 
the eq. (3). Our task is to find out how accurately the lens parameters can be determined. 
Instead of a massive Monte-Carlo simulation we adopt another approach. First, the values 
of Fq, Aq, fs, Umin, imax, ^0 are fixed. The first two are assumed to be known to the 
observer, while the latter four are to be determined. The x^ of the fit between a model and 
a string of quasi-data points {Fdata,k,tk), k = 1, 2, 3, ..., n is approximated with the sum: 

X2 V"^ (,-f model ,k ^ data,k) ~r ^^ /_ n 

- Z^ XI ' y^'^) 

k=l ^k 



model, k 



{Fq - F,o) + F,,fc, (56) 



where Fmodei,k are the intensities which follow from the model we are fitting, the 
quasi-observations Fdata,k are given with the eq. (2), and the measurement errors A^ 
are given with the eq. (4). Our task is to determine the values of four parameters: 
cti,2,3,4 = {fs, Umin, ^max, ^o) by minimizing X^. Note that because of measurement 
errors the value of X^ can never be zero. It is minimized for the correct choice of the 
four parameters, when its value is n. We shall also determine the confidence ranges of 
the four parameters for which the value of X^ is within a chosen range of its minimum. 
For simplicity we assume that n equally spaced measurements cover the time interval 

For small errors we can expand X"^ around its minimum to obtain: 

4 4 

A^' -EE A,, Aa,Aa„ (6) 

7=1 i=l 
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where Actj = amodei,i — ctdata,i ii = 1, 2, 3, 4) is a vector of the parameter differences between 
the fitting and the model values, and 
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at ftj = adata,i7 {i = 1? 2, 3, 4). For a quadratic expansion the confidence ellipsoid it describes 
scales linearly with the measurement errors. The errors with which all model parameters 
are determined scale linearly with the quantity a = Aon~^/^, for n ^ 1 and Aq ^ 1. 

The best MACHO and OGLE microlensing events had n ~ 60 and Aq ~ 0.04, for 
the corresponding a ^ 0.005 (cf. Udalski et al. 1994a, Alcock et al. 1995). A more 
typical values were (n, Aq, a) ~ (30,0.1,0.02). However, with the introduction of effective 
follow-up observations by the PLANET (Albrow et al. 1996) and OMAN (Pratt et al. 
1996). a substantial increase of n, a reduction of Aq, and the corresponding reduction of 
the parameter a have already been achieved in some cases, and farther improvement is 
expected in the near future. 

In the following section we investigate some problems with the determination of model 
parameters imposed by the finite value of the a parameter. 



3. The Degeneracy 

There are two regions in the parameter space for which there is a near degeneracy, i.e. 
with any realistic accuracy of the measurements it is not possible to determine the unique 
values of all model parameters. 

The first troublesome case is that of a large impact parameter. In the limit Umin ^ 1 
the eqs. (2) and (3) may be transforemed as follows: 



A 



2 

1 + ^ 



^2 
"mm 



U- \ut,r, + (t/to)' 



|2' 



for M ^ 1, tr 



0. 



Within this approximation we can write (cf. eqs. 2-4) 
Fit) ..„.„/. 2 



:i - /.) + /s X 1 



2/s 



-u: 



+ (t/to) 



NLn + {^lUff 



It is straightforward to verify that substituting the parameters 



/s,l = JsC 



U 



miriA 



^rn.iri\-^ , f^i 



0,1 



tnC 



-1 



(9) 



(10) 
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into the eq. (9) we recover the same formula for the Fs/Fq. Therefore, if the set of 
parameters: (/s, Umin, to) is a solution then the set {fsC*, UminC, toC^^) is also a solution, 
where C is an arbitrary positive constant which satisfies fgC^ < 1. 

Analytical considerations seem to imply that the trouble sets in when the impact 
parameter is very large, Umin ^ 1- In reality the situation is much worse. An example of 
the problem is shown in Fig. la, in which the two light curves shown with a solid and 
a dashed line, respectively, are almost identical. The truly troublesome aspect of this 
example is the value of the impact parameter: Umin = 0.5. This implies that in practice the 
degeneracy covers a broad range of impact parameters, and a seemingly robust microlensing 
event with the peak magnification in excess of 2 may have a blend which is photometrically 
undetectable in the event's light curve. 

Another troublesome case is when the minimum flux, i.e. the fiux measured far from 
the microlensing event, is dominated by a blend or by any background. This is a generic 
case in very crowded fields, like the bulge of M31, in which microlensing events are to be 
detected by the "image subtraction" method (Bouquet et al. 1996, Crotts & Tomaney 1996, 
Tomaney & Crotts 1996, 1997). As the lensed object contributes little to the total light at 
minimum, the microlensing event can be detected only when the peak magnification is very 
large, i.e. the impact parameter is very small. In this case we have Fsq/Fq <^ 1, Amax ^ 1, 
Umin *^ 1, and a significant change in the brightness occurs only close to tmax when the 
projected distance between the source and the lens is very small, i.e. when -u <^ 1. We have 



A« 1 + ->1, 

u 


for u < 


1 1 •' ^ — ^ 1 


fs 


u Uj^ . 

."'mm 


+ it/toy]'^' 



u 

and therefore 

Fit) 

-±l = {l-f,) + f,A^ 1 + ^ = 1 + — ^^— :;tt7^, for /^ « 1, tmax = 0. (12) 

-To 

It is straightforward to verify that substituting the parameters 

Js,l = Js^j Umin,! = Umin^ i ^0,1 = ^0^ ) (l<j) 

into the eq. (12) we recover the same formula for the Fs/Fq. Therefore, if the set of 
parameters: {fs, Umin, to) is a solution then the set (/^C, UminC, toC~^) is also a solution, 
where C is an arbitrary positive constant which satisfies fgC < 1. An example of such a 
case is shown in Fig. lb. 

A more quantitative way to analyze the degeneracy is to plot the confidence contours 
in the to — Umin and in the fg — Umm planes. Two sets of such plots are shown in Fig. 2; 
they correspond to the two cases presented in Fig. 1. It is apparent that the contours are 



-6- 



highly elongated, i.e. a certain combination of model parameters may be determined with 
a reasonable accuracy, while the other combination is subject to a very large error. The 
photometric accuracy parameter a = Aq/v^ = 0.007 was assumed, which may be achieved 
for example when Aq = 0.05 and n = 50. 

The time of maximum magnification, tmax, and the peak brightness, Fmax are not 
subject to any degeneracy in their determination as long as there are many points within 
FWHM of the light curve. Unfortunately, these two parameters are of no particular 
significance for any inferences from the microlensing observations. Those parameters which 
are significant: the event time scale tg and the impact parameter Umin tend to have strongly 
correlated errors, as shown in Fig. 2. It is a general property of blended microlensing events 
that some combination of the (to? ""mm) values can be determined with a good accuracy, 
while some other can be determined only poorly. In both cases of degeneracy as described 
with the eqs. (10) and (13) the product to >< Umin can be measured well, but the ratio 
to/umin cauuot bc determined from the observations, as it can be equal to any number. 
Therefore, we cannot determine the values of to and Umin in these cases, only the value of 
the product of the two parameters. Note, that in accordance with this analytical reasoning 
the confidence countours as shown in the to — Umin plane in Fig. 2 appear to be streached 
along a hyperbolic shape. 

There are two dimensionless parameters which determine the shape of a blended 
microlening light curve: the impact parameter Umin, and the fraction of minimum light 
contributed by the lensed source, fs = Fso/Fq. For a chosen value of the accuracy parameter 
a = 0.007 we calculated values of the standard deviations in the determination of all four 
parameters: F^o, Umin, to, tmax as a function of Fsq/Fq and Umin, and they are presented 
in Fig. 3. The lines are labeled with the values of standard deviation. For example, in the 
upper left corner the error in the determination of the fraction of light in the lensed object, 
SFso/Fso is equal 1.0 (i.e. a huge error) along the uppermost line, and the error is modest 
at 0.1 along the lowest line. If the accuracy parameter is smaller by a factor 10, i.e. if 
a = 0.0007, then the labels along all lines would be reduced by a factor 10, i.e. the accuracy 
of the determination of Fsq/Fq would be 0.1 along the uppermost line, and 0.01 along the 
lowest line. A reduction of the a parameter by a factor 10 can be accomplished by reducing 
the photometric error Aq by a factor 10, or by increasing the number of photometric 
measurements n by a factor 100. 

An inspection of the Fig. 3 reveals that only the time of maximum magnification t^ax 
can be well determined over a significant fraction of the parameter space. The other three 
parameters: F50, «„«„, to are virtually impossible to measure above the lines with the label 
1.0. Note that our choice of the accuracy parameter a was optimistic by the standards of 
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current microlensing seraches. 

As the confidence contours are so elongated (cf. Fig. 2) it is interesting to check if 
the best combination of the interesting parameters can be measured with a significantly 
higher accuracy than possible for each of the parameters separately. This is shown in Fig. 
4, where the right two panels present the lines labelled with the corresponding values of 
a standard deviation for the best linear combination of {umm, F'sq/Fq) (upper right panel) 
and {umin,to) (lower right panel). Note, lines in these panels are placed much higher than 
the corresponding lines in Fig. 3, i.e. the errors are strongly reduced by choosing the 
best combination of the parameter pairs, corresponding to the short axis of the elliptical 
confidence contours in Fig. 2. On the other hand, the opposite choice, corresponding to the 
long axis of the elliptical confidence contour, and presented in the two left hand panels in 
Fig. 4, reveal errors even larger than those corresponding to each parameter separately, as 
presented in Fig. 3. 

Let us now ignore the presence of blends, and let us follow the standard procedure, 
assuming that Fso = Fq. With this assumption we determine the values of the event time 
scale tg and the impact parameter m^j^. The parameters so determined may be compared 
with their true values to and Umin in Fig. 5. It is clear that by neglecting the blend we 
always overestimate Umin and we underestimate to, with the error increasing towards lower 
values of Fsq/Fq, i.e. stronger blending, and lower values of impact parameter. 

With the blending effects being so important we should check when the presence of 
the blend can be established by means of photometry alone. It was done so with the two 
double lenses (OGLE #7, Udalski et al. 1994b; DUO #2, Alard et al. 1995) and with one 
single lens (OGLE ^5, S. Mao, private communication 1995, Alard 1996b). Unfortunately, 
this task is difficult for events caused by a single point mass, as shown in Fig. 6, where the 
change in the value of X^ is shown as a function of two dimensionless parameters: Fsq/Fq 
and Umin- This is the difference between the X^ values of a blended event as fitted with 
and without a blend. If the change is large then photometry alone can clearly demonstrate 
that the lensed source has a blend. The parameter space above the thick line AX^ = 2.8 
is photometrically degenerate: in this region it is not possible to detect the presence of a 
blend with the photometric accuracy parameter a = 0.007 at the confidence level 90%. 



4. Discussion 

The main result of this paper is somewhat discouraging: with the currently typical 
photometric coverage and accuracy the light curve of a single microlensing event cannot be 



used to determine the presence of a blend unless the impact parameter is small, Umin < 0.3, 
and the blend may be photometrically undetectable even if the impact parameter is very 
small (cf. Fig. 6). In the real world we do not know the precise value of the measurement 
errors, and there is a significant contribution to the errors which is non-gaussian (Udalski 
et al. 1994a). Therefore, the photometric border of blend detectability as indicated by 
the thick line in Fig. 6 is optimistic, with the true border located at even lower values of 
Umin- Yet, the blending must be very common as demonstrated by the two double lensing 
events (OGLE #7, Udalski et al. 1994a; DUO #2, Alard et al. 1995), and by the recent 
analysis of the DUO results (Alard 1996b). Unless the blending is somehow corrected 
for, the microlensing time scales to are underestimated, and the impact parameters are 
overestimated (cf. Fig. 5). This implies that the lens masses and the optical depth are 
systematically underestimated. 

There are various ways in which the blending may be detected or at least statistically 
corrected for, as discussed by Alard (1996b). In case of a double lens, with caustic crossings 
indicated by the light curve, a blend reduces the apparent magnification between caustic 
crossings below the theoretical minimum value Amin = 3 (Witt & Mao 1995) - this was 
apparent in OGLE ^7 and in DUO ^2. In case of any event, double or single, a blend can 
be detected through a correlation between the image centroid and the apparent brightness, 
as first noticed in DUO ^^2 (Alard et al. 1995). This is possible if the angular separation 
between the lensed object and the blend is not too small, presumably no less than half a 
pixel, or so. If the blend has a different color than the lensed star then it may be uncovered 
with multi band lightcurves of a microlensing event (Buchalter & Kamionkowski 1996). If 
the luminosity function well below the detection threshold is known, say from the HST 
imaging, then the blending can be statistically corrected for by placing artificial stellar 
images on the CCD frames. 

The last approach has its limitations. As long as we do not know the distribution 
of separations and luminosity ratios of binary stars in the target population we can only 
put artificial stars randomly on a CCD frame. But this procedure will underestimate the 
number of very close binary pairs. Fortunately, if the detected number of microlensing 
events is very large we can use the microlensing itself to determine the statistics of binary 
stars, as they will produce repeating microlensing events (Di Stefano & Mao 1996). 

This work was supported with the NSF grants AST-93 13620 and AST-9530478. 
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Fig. 1. — Two examples of degeneracy are shown, (a) A large impact parameter. The 
unblended event with Umm = 0.5 and to = 1.0 (solid) is well fitted by the model with 
'^'min = 0.4, tg = 1.15 and Fsq/Fq = 0.73 (dashed). The difference in the goodness of fit 
as defined with eqs. (5) is AX^ = 0.52. (b) A strong blending. The solid line shows the 
microlensing light curve with the parameters: Umin = 0.15, to = 1-0, Fsq/Fq = 0.07. The 
dashed light curve is the best light curve for unblended event with -u^j^ = 0.97, tg = 0.30. 
The difference between the two corresponds to AX^ = 2.63. In both plots a large dot with 
the error bar indicates Aq, the photometric error at minimum light. 
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Fig. 2. — Confidence contours in the (to — Umin) and {Fsq/Fq — Umin) planes for two cases of 
degeneracy from Fig. 1: a large impact parameter {Umin = 0.5, upper) and a strong blending 
{Fso/Fq = 0.07, lower). Confidence levels are 68% and 90% for one parameter of interest. 
The photometric accuracy parameter a = Aq/v^ = 0.007 was adopted. 
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Fig. 3. — The accuracy of the best fit parameters (68% confidence) as a function of Umin 
and Fsq/Fq. The hues correspond to a constant error in the best fit value determination for 
a given parameter. For Fsq/Fq and to the fractional errors are shown. The errors in tmax 
are in units of to, while the errors in Umin are in units of Einstein radius. The lines are 
equally spaced in error value, between 0.1 to 1.0 for Fgo, Umin and to, and from 0.01 to 0.1 
for tmax- All parameters are best measured when the impact parameter Umin and the blend 
contribution {Fq — Fsq)/Fo are both small, i.e. in the lower right hand corners of every panel. 
The photometric accuracy parameter a = Ao/y^ = 0.007 was adopted. The ploted errors 
scale linearly with a provided that the number of photometric measurements n is large and 
the photometric error Ao is small. 
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Fig. 4. — The same as Fig. 3 but for two linear combinations of {umm, -^soZ-^o) (upper 
panels) and {Umin,to) (lower panels), with the two combinations corresponding to the long 
axis of confidence ellipse (the largest errors, left panels), and the short axis of the ellipse 
(the smallest errors, right panels). 
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Fig. 5. — The effects of fitting a standard light curve (without a blend) to the blended 
microlensing events. The plots show the ratio of the best fit value of m^j^ (upper panel) 
and tg (lower panel) as a function of \og{Fso/Fo) and Umin- If blending is ignored then the 
impact parameter Umin is overestimated and the timescale to is underestimated. The lines 
are labeled with the values of u'^^^/umin and tg/to, and are equally spaced in these values. 
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Fig. 6. — The difference AX^ between the best fit model with F^o = Fq and the best fit 
model with adjustable Fgo is shown as a function of the impact parameter Umin and the 
fraction of the lensed light Fsq/Fq. Solid lines correspond to AX^ = 0.01,0.1, 1.0, and 10.0 
for the photometric accuracy parameter a = Ao/y^ = 0.007. It is impossible to discriminate 
statistically between the two types of models for AX^ < 2.8, i.e. above the thick line (90% 
confidence) . The two areas of strong analytical degeneracy, one for large Umm and the other 
for small Fsq/Fq, are two parts of a single very large region. When Fsq/Fq approaches 1 the 
weak blend cannot be detected from a light curve alone, but in this case the blend is of no 
practical consequence. 



